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Abstract 

The introduction of an explicit negative radial excitation contribution in the 
hadronic side of the light cone QCD sum rule (LCSR) of Belyaev, Braun, Khod- 
jamirian and Riickl, can explain the large experimental value of gr>* Dm recently mea- 
sured by CLEO. At the same time, it considerably improves the stability of the sum 
rule when varying the Borel parameter M 2 . 

PACS: 14.40. Lb (Charmed mesons), 13.75. Lb (Meson- meson interactions), 11.55.Hx (Sum rules) 



1 Recalling the difficulty 

The discrepancy of the LCSR prediction for gD*Dn with the recent experimental result by 
CLEO is a striking and puzzling one. We recall the reader that the LCSR prediction made 
in ref. [ij reads gD*Dir = 12.5, which got even smaller after the radiative corrections have 
been included 0, namely gD*Dn — 10.5 ± 3.0. On the other hand the experimental value 
is gD'D-K = 17.9 ± 0.3 ± 1.9 0. Thus the discrepancy is large 1 . The above LCSR result 

* Unite propre de Recherche du CNRS - UPR 7061 
tUnite mixte de Recherche du CNRS - UMR 8627. 

1 Notice that we speak here of the discrepancy of central values. Of course, there are errors in the sum 
rule estimate, but it is difficult to draw strict conclusions from them, since their estimate is itself, of course, 
uncertain. 



has been very carefully discussed, the calculation has been improved several times, and 
it has also been verified by using other sum rule techniques [|J. Although the QCD sum 
rule approach certainly suffers from large uncertainties, in most cases a good agreement 
with experiment has been obtained. Therefore, we cannot be satisfied with concluding 
by a comfortable attitude of skepticism against the whole sum rule approach. Still more 
puzzling is the fact that this quantity does not seem a priori to have anything particularly 
exotic. In addition, the other theoretical approaches do not show such a discrepancy for 
this quantity. A careful discussion meant to reduce the uncertainties presented by quark 
models, and performed in the framework of Dirac equation ||, prior to the experimental 
measurement, has led to a result gD*Dn — 18. It should be stressed that this result has 
been obtained in the heavy quark limit (for the latest quark model discussion of go*Dn see 
ref. !|). The recent (quenched) lattice QCD calculation has led to go'D-n = 18.8±2.3l^o [fl- 
it is therefore important to understand the specific difficulty which the standard sum rule 
approach seems to encounter in this case 0. In ref. || it has been noted that the simple 
quark-hadron duality ansatz which works in the one-variable dispersion relations might be 
too crude for the double dispersion relation. 



2 An appealing solution: a large negative radial exci- 
tation contribution to the LCSR 

In sum rule calculations it is usual to assume that the higher state contributions can be 
included in the perturbative estimate of a continuum 2 . In other words, one does not include 
an explicit contribution of an isolated excited state. The first reason is probably pragmatical: 
it seems unnecessary to recourse to large excitation contributions if the stability criterion 
can be well satisfied without them. The other, more theoretical, reason is that the excitation 
contributions are exponentially suppressed with respect to the lowest state by the Borel 
procedure, so that one can sufficiently account for them by the rough procedure of the 
perturbative treatment of the continuum. Indeed, in the LCSR approach of ref. 0, this 
point of view has been adopted. We suggest, instead, that this neglect of an explicit 
radial excitation contribution may be the origin of the above discrepancy between the 
LCSR prediction and the experimental value for go*Dn- Of course, we are aware that, 
as explained in Belyaev et al. ||, in their method such contributions are suppressed by 
the Borel exponential. Our claim is that, in this particular case, the Borel suppression 
is not sufficient to allow to neglect them. The inclusion of the explicit radial excitation 
contribution to the hadronic side of the LCSR (often referred to as the left hand side - 
l.h.s.- of the sum rule) also offers an appealing explanation for the failure of the LCSR 
prediction. 

A first indication in favor of this new proposal is that, after including a radial excitation, 
the stability of the sum rule, under the variation of the Borel parameter M 2 , is improved. 
This criterion is however subject to an uncertainty since we do not really know what is the 
accuracy of the calculation of the theoretical r.h.s. It would be therefore good to present 
cross-checks of our assumption, which is what we will do in the next sections. 

2 Except for the critical case of quark masses in the pseudoscalar correlator method, or in general for 
some refined calculations like the ones discussed in ref. [fioi . 
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Before embarking on those issues, let us first explain why the assumption of a large 
radial excitation coupling to the sum rule, with a negative sign relative to the ground state 
contribution, seems so interesting. Since the purpose of this letter is to communicate our 
proposal, and not to make the best numerics, we leave apart for the moment the radiative 
corrections and write the sum rule of Belyaev et al. as: 

_ f(M 2 ) 

gD*Dir — -~r~ f — , l-U 

JDJD* 

where the function f(M 2 ) is the r.h.s. of eq. (44) of ref. 0, namely 

^ - hum M ' 2 exp i^ 3 ^) Y mVh " - e-,o/ i + • • • (2) 

We do not write the higher twist terms since it would make the expression lengthy and 
would not help understanding our proposal. Those terms are numerically very important 
and are indeed included in our calculation. It is then found that, within the Borel window 
2 GeV 2 < M 2 < 4 GeV 2 determined by the standard criteria, the function f(M 2 ) is 
monotonously decreasing, and the variation is as large as 20%. Therefore, there is no truly 
good plateau. The authors quote an average 

J( M % GeV*<M 2 <4 GeV* = ' 51 ± - 05 G ^ > ( 3 ) 

which then yields the central value : 

9D'Dn = 12.5 , (4) 

indeed much too low. Let us now introduce a radial excitation contribution to the hadronic 
l.h.s. of eq. (44) of ref . 0, or equivalently write: 



gD*Dn 



D* 



m 2 D , — m 2 D \ ( m 2 D ,, — m 2 D * 



f(M 2 ) - Rjy exp ^7777^ " R »*' ex P 



2M 2 J 1 \ 2M 2 



(5) 



Note that we have two extra contributions: either the D or the D* is excited; they are 
denoted as (D',D*) and (D,D*'), leading respectively to : 

/ m D , \ 2 m D *, 

tiD' — JD'jD*gD*D'-K , tlD*> — J D J D*' g D* 1 Dtt ■ (0) 

\m D J m D * 

We assume that the higher (D',D*'), (D",D*), ...contributions are still included in the 
continuum part of the model. The completion of the procedure requires also a new value 
of So, since the lowest radial contributions are no more comprised in the continuum part, 
but are handled separately. 

For simplicity, since there is no sense in requiring too much precision, we will assume: 

m D , = m D *, , f D , = fu*i , (7) 

the spin-spin effect being expected to decrease for higher states. We also assume that 
gD*D'ir = gD*'Dn = g', which is more questionable, since the ground state is present, and D 
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differ appreciably from D*. This, however, should not alter the qualitative conclusion that 
we obtain. Note also that the difference between the 0~ and 1~ states would disappear in 
the heavy quark limit. With these assumptions, we have 



Rd> m D *im D * f D * 



Rd*> rn D f D 



(8) 



This equation (§) relates the two radial contributions since Rd'/Rd* 1 is directly calculable 
if we know the mass tub* 1 - The exponentials accompanying Rd'/Rd*' are the trace of the 
Borel suppression of radial excitations; they are increasing functions of M 2 . Therefore, 
if we choose a negative Rd',d*', we can compensate for the decrease of f(M 2 ) and thus 
simultaneously: (i) improve the stability of the sum rule and, (ii) increase the magnitude 
of gD*Dn (see fig. [I] for illustration). We have of course to take a rather large radial 
contribution in order to get a significant effect. To avoid adding more freedom, we fix the 
excitation spectrum by a model calculation which is completely independent. Using Dirac 
equation, with numbers taken from our previous treatment ||, we find in the heavy quark 
limit that the excitation energy is 0.5 GeV, and therefore the mass of the excitation is 
2.5 GeV. Let us recall that the large width ||TT| expected for such a state makes it difficult 



to observe. 

The remarkable fact is the following: taking a negative Rd*i (Rd 1 ), and imposing that 
the variation of go'Dw is 5% or less in the allowed window for M 2 , (instead of the previous 
20%, in the same Borel window), we find that Rr>*' must be in the range 

- 0.25 GeV 2 < R D *, < -0.1 GeV 2 . (9) 

Then, depending on the value of Rd*> in the above range (|9]) the predicted value of gn*Dn 
varies as 

17 < g D ^ < 25. (10) 

It increases with |-Rd*'|. It is much larger than before. We obtain at the same time a much 
better stability in M 2 , and a much larger coupling constant, just in the desired range to 
agree with experiment. 

For Rd*i = —0.15 GeV 2 , we would obtain a perfect stability in M 2 (accuracy of the 
order of 1%) and the value go'Dw = 19.4, that is the value found in lattice QCD. But at 
the present stage, there is no sense in trying to get such an accuracy in stability, given the 
very rough precision of the calculation. 

For the present note, we have adopted the parameters of ref. [0 for the allowed range 
of variation M 2 ; we choose also to keep the continuum threshold to be the same as sq, the 
continuum threshold of the two-points annihilation sum rules proposed by these authors, 
although this may be questioned since we now take apart certain contributions of excited 
states; we do not indulge in any fine tuning. We have also not allowed the annihilation 
constants fo, fD* to vary with M 2 , but instead, we have taken the final value given in 
the same reference ]IJ: indeed, we should check whether fo ,fo* as well as gD*Dn are 
independent of M 2 . As we stated above, the inclusion of the radial excitation stabilizes the 
LCSR. This does not make the same effect in the sum rules for fjj, f^,*. As we shall see, the 
radial excitation contribution is neither stabilizing the sum rule (w.r.t. the variation of the 
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Figure 1: The effect of inclusion of the radially excited .OW'-meson contributions to the hadronic part of 
the LCSR: In the lower figure the result of ref. is reproduced (no explicit radial excitations considered) ; 
In the upper curve we include the radial excitations as indicated in eq. (|^) with the value of Rd*i varied in 
the range specified in ([)]) . The latter range arises from the requirement of the stability of the sum rule (0) 
to be within the 5% level, when the Borel parameter is varied as 2 GeV 2 < M 2 < 4 GeV 2 . 
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Borel parameter M 2 ), nor it creates a difficulty: it simply gives a negligible contribution, 
and does not modify the previously estimated values for fn, fr>*. 

We note that the minus sign for the radial excitation contribution is crucial, and we 
would like therefore to have some argument for this from models. There is no very trustable 
model for radial excitation strong decay at the quantitative level. Nevertheless, it is in- 
dicative for our purpose that the rather standard model for strong interaction couplings, 
the non relativistic quark pair creation (QPC) model, used for calculational simplicity 
with harmonic oscillator wave functions, give precisely a stable negative sign for the prod- 
uct f d f d*' g d*' Dtt relative to the similar product for the ground state fnfD*gD*Dn- Hence, 
since, in the sum rule, the latter product is positive, this means that Rd*>, which is precisely 
fDfD*'gD*'Dn up to positive mass ratios, should be negative, as we have found in the sum 
rule. 



3 Check of D —> n£u semileptonic decay 

We note that exactly the same quantity Rr>*> will appear in a t-channel analysis of D — ► m! 



decay. Indeed, using an unsubtracted dispersion relation, as required for the form factor |L2 
the addition of the residue of the pole corresponding to the radial excitation results in : 

p C 2 ) = 1 fp*9D*D7T 1 m D *R D *> 

+ 2m D * 1 - q 2 /m 2 D * 2f D m 2 D „ 1 - q 2 /m 2 D „ 

So the question is whether our "large" contribution of radial states in the LCSR leads to 
effects in the D — > niu which are compatible with the facts. The answer is first that its 
impact is rather small, and anyway rather in the right direction. Indeed, we find that for 
Rd* 1 — —0.15 GeV 2 , allowing go'D-n — 19, the residue of the radial excitation is almost one 
order of magnitude smaller than the ground state one : 

F+iq) ~ 1 - 0.25g 2 " 1 - 0.16g 2 + - (12) 

It should be noted that the ground state contribution alone would lead to a too large 
rate (the numerator should be about 0.8 at most instead of 1.15 to get the correct rate 
with only the first pole). Therefore, the negative radial contribution is not worrisome 
but, quite the contrary, it is not sufficiently large to balance the excess of the lowest state 
contribution. Other radial excitations, as well as the continuum, must contribute, which is 



in agreement with the conclusions of ref . . Analogously, a negative contribution appears 



in the Becirevic-Kaidalov model [13], from remote singularities. 



4 Check of annihilation constants and Adler— Weisberger 
sum rule 

Another worry could be that this large contribution could contradict other sum rules con- 
cerning either the annihilation constants or the strong couplings. The annihilation constants 
fo and fo* are determined from the QCD sum rules by using the P-P and V-V correlation 
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functions, respectively. At the same time, an Adler-Weisberger sum rule constrains the 
strong interaction coupling constants. Since the latter is stronger, as we shall see, we begin 
with it and write it, in the case ofn + D^n + D scattering, as : 



J IT 



n-0 '"'/V ' 



9D*( n )Dw\ + orbitally excited states = 1, (13) 



where we have made explicit only the D* ground state contribution and the one of its radial 
excitations ; the allowed orbital excitations are the + and 2 + . In ref. f|, we have estimated 
these contributions of the ground state and lowest orbitally excited states (D**) to be 70%, 
whence 

< Vox = 0.55 . (14) 

Now, since we have \R D (*),\ > 0.1 GeV 2 , we can deduce a lower bound on fa*/ = f^i from 
eq. (|) 

\fjy\> 0.02 GeV. (15) 

This is one magnitude smaller than the ground state annihilation constant. Studying the 
QCD sum rules for fp, fc*, we find that the resulting contribution from the radially excited 
states is very small at the level of this lower bound. For the pseudoscalar sum rule, one has 
to change the l.h.s. according to 

,2 ,2 . f m D'\\ 2 ( m D> - 

Numerically, that change is completely negligible: it adds less than 0.0016 GeV 2 (times 
the exponential factor, dependent on M 2 , and which is smaller than 0.25) to /|> = 0.17 2 = 
0.029 GeV 2 . Therefore, we have much room left even if the bound (|T4]) were to be lowered. 
For the vector sum rule, the conclusion is quite similar. The substitution to be made is 

SI- - SL + feV SI; exp (-2^2^) . (17) 



m 



D 



The change is numerically even smaller in this case. 

In contrast to the LCSR for fDfD*9D*Dir, the radial excitation term does not improve 
the stability of the sum rules for the annihilation constants: indeed, the r.h.s. is decreasing 
with M 2 , and the contribution from the radially excited state increases, since it is positive. 
But it does not deteriorate the stability, since it is very small. 



5 Conclusion 

The result of the proposed modification of the sum rule calculation is very encouraging. 
Introducing the radial excitation contribution significantly improves the value of gD*Dn, 
and at the same time the stability of the sum rule with respect to the Borel parameter 
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M 2 . This means that the effect of such a radial state is not properly accounted by the 
standard perturbative continuum contribution. Unless one imposes an unreasonable degree 
of stability, one is not able to fix the magnitude of the radial contribution accurately. The 
latter may vary by more than a factor of two, which induces an important variation of 
9D*Dn, around 30%. This should not be viewed as a problem since our main goal is to prove 
that a large value for gn*Dn, i-e. in the range allowed by experiment, is possible in the 
LCSR approach. At the very least, we can say that the old solution with the low value of 
9d* Dn and no explicit radial excitation, is not compelling, and alternative ones with large 
values of gD*Dn are favored. Of course, one should wonder whether the success survives 
when the calculation of the theoretical side is improved. In particular, does it survive the 
introduction of the radiative corrections, which are large? Our answer is positive, but we 
would like to present it within a more extensive discussion of the parameters involved in 
the sum rule calculation, especially because the threshold parameter s should depend on 
the fact that we separate an explicit contribution of some excited states. 

The main reason for having such a large effect of radial excitations without exceedingly 
large couplings is that the Borel exponential suppression effect is very weak in the LCSR: 
at most a factor 0.5 (at the lower end of the allowed range) for the relative suppression of 
the radial excitation with respect to the ground state. The Borel suppression is much less 
effective than in the standard sum rules, e.g. for correlators leading to /dm. This is due to 
the additional factor of 1/2 in the exponent and to the fact that M 2 is twice larger on the 
mean in the allowed range (2 GeV 2 < M 2 < 4 GeV 2 , instead of 1 GeV 2 < M 2 < 2 GeV 2 ). 
At the lower end of the range, the suppression is around eight times less effective in the 
LCSR. 

Of course, it would be good to examine whether one gets similar improvement for the 
other sum rule approaches to gn*Dw On the other hand, the situation encountered here 
underlines the usefulness of an alternative to the Borel transformation method, able to 
isolate with more efficiency various states. In some cases, the FESR proposed for the case 
of the quark mass sum rules by Kambor and Maltman |14| could be a good alternative. 
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